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Abstract 

In this paper we establish global existence and uniqueness of the 
solution to the three-dimensional Vlasov-Poisson system in presence 
of point charges in case of repulsive interaction. The present analysis 
extends an analogeous two-dimensional result pQ. 



1 Introduction 

In this paper we study the time evolution of a three-dimensional system con- 
stituted by a continuous distribution of electric charges, a plasma, coupled 
with N charged point particles. 

All the charges, as well as the plasma, have the same sign so that the 
interaction is repulsive. For simplicity we assume that the charges and 
the masses of the point particles are unitary. If / = f(x,v,t) denotes the 
mass distribution of the plasma and {£ a }a=i are the positions of the point 
particles, the dynamics of the system is described by the following system 
of equations 



dtf + v- V x f + (E + F) ■ V v f = 



E{x,t) 



p(x,t) 



y 



r3 F - y\ 



P(y,t)dy 



N 



f(x,v,t)dv 
x - £ a (t) 



L{t)=Va(t), a = l,...,N 



(1.1) 



V 
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In absence of the charges, the system (jl.ip reduces to the well-known Vlasov- 
Poisson equation, which has been widely investigated in the last years. 

The difficulty of the Cauchy problem associated to the Vlasov-Poisson 
problem increases with the dimension of the physical space. 

In two dimensions, satisfactory existence and uniqueness results go back 
to [8] and [5]. The three-dimensional problem was solved in the nineties in 
[2] j PS] j [II] > by a careful analysis of characteristics associated to the Vlasov- 
Poisson system (Lagrangian point of view) or by estimating the moments of 
/ via a more genuine PDE technique [7] (eulerian point of view). We address 
the reader to the monograph [2 J for a complete analysis of this equation and 
additional references. 

When point charges enter in the game, the situation changes drastically 
even if complete repulsivity is assumed. The extra singular field could, in 
principle, produce extremely large velocities of the plasma particles and, in 
turn, a large spatial density and a blow-up of the electric field in finite time. 
However we know that it is not the case in dimension two. Indeed, in pQ it 
was shown a global existence and uniqueness result for the two-dimensional 
Cauchy problem associated to system (jl.ip . The basic ingredient of the 
analysis in [1] was the introduction of the energy of a trajectory of the 
plasma in the reference frame of a suitable point charge. The advantage of 
using this energy function is two-fold. From one side it controls the motion. 
From the other, its time derivative along the trajectory does not involve the 
singular part of the electric field. Combining this idea with the well-known 
fact that, in dimension two, the electric field generated by the plasma is 
linearly bounded by the maximal velocity of the particles, one can then 
conclude (see [1] for the details). 

Our purpose here is to study the more complex three-dimensional prob- 
lem. Our approach relies heavily on the adaptation of the method in [9], |10j 
and [11] to the present situation, for which the energy associated to a tra- 
jectory (see (|3.7p or (|4.3p ) plays an essential role. We explain the main steps 
of this adaptation in the next section after some preliminary considerations 
have been presented. 

We mention that a related issue concerning the Cauchy problem for 
the Vlasov-Poisson equation, namely the instability arising from a singular 
perturbation of the field, has been recently analyzed in [6], see also [3] or 
[3]. In that situation the extra singular field is due to reflecting boundary 
conditions. 

The plane of the paper is the following. In Section [3] we solve the Cauchy 
problem for one single charge (N = 1). This result is well suited to be easily 
extended to the full problem N > 1. This extension is done in Section HJ 
Section [5] is finally devoted to comments and criticism. 
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2 Preliminaries and general strategy 

We use this section to fix the notations, to recall some well known estimates 
(see e.g. [11], [2] and references quoted therein) and to illustrate the strategy 
we follow to treat the present problem. 

Let / = f(x, v), (x, v) € M 3 x M 3 , be a probability density (namely / > 
and f f dxdv = 1). We denote by 



the spatial density. 

Also, K and Ki,i = 1, . . . , will stand for positive constants depending 
only on the kinetic energy ()2.ip and on , which we assume to be finite. 

Moreover C will denote any numerical positive constant. 

For any M > 0, we have 




(2.1) 



the kinetic energy and by 





Optimizing in M, we find 




3/5 



whence 



p\\ L S/3 < K ± . 



(2.2) 



Moreover, defining for R > 




we have by using Holder inequality 



/ 



Pr{x') 



dx' = 




x — X' 



./|2 
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Optimizing in M and using (|2.2p we obtain 

%dx' <K 2 \\p R fll. (2.3) 



\x — x \ 

Summarizing the previous considerations we are led to the 

Proposition 1. There exists a constant K > (depending only on J \v\ 2 fdxdv 
and || /|| l°° ) for which 

-PK^l-dx' <KR^ (2.4) 
\x — x \ z 

for all R>0. 

Proof. Estimate (|2.4[) follows from (|2.3p . realizing that ||pi?||oo 

□ 

We now come to Problem (jl.ip and define the class of solutions we wish 
to deal with. 

Let {(£,aO,Vao)}a=i denote the initial positions and velocities of the point 
charges. Let /o be a compactly supported probability density satisfying, for 
some positive 5o, 

min{|x - £ q0 | I (x, v) G supp(/ ), a = 1, . . . , N} > S . (2.5) 

Let T > 0. We say that • • • , £at, /) is a solution of Problem (jl.ip on 
[0, T] with initial datum (£ Q o> VaO, a = 1, . . . , N, /o) if 

^GC 2 ([0,T]), /GT^n^nr), p G (|0, T], L°°) (2.6) 

and for i G [0, T] we have 

'£*(*) =??a(i) 



Moreover 

f (X(x,v,0,t),V(x,v,0,t),t) = f (x,v), (2.8) 
where for all r,t G [0, T] 

(X,V)( V! r,t) :R 3 \U a {Ur)} xR 3 ^R 3 \U a {Ut)} xl 3 

is an invertible flow such that 

\X(x, v, r, t) - Ut)\ > S(T), V(x, t;) G supp(/(r)) (2.9) 
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for some S(T) > 0. It satisfies for t G [0, T] 
( d 

— X(x, V, T, t) = V(x, V, T, t) 

at 

d \ - X(x,V,T,t) -£ Q (t) 

-V(x, v, r,t) = E (X(x, v, r, t), t) + g \ x{x ^ t) _ UtW 

(X, V) (x, v, t, t) = (x, v) G supp(/(r)). 

(2.10) 

Given p G L°° ([0, T], L 1 n L°°(IR 3 )) it is a well-known fact (see e.g. [7]) 
that the corresponding field E belongs to L°° ([0, T] x M 3 ). Moreover, it is 
almost-Lipschitz in the sense that for all (x, y, t) G M 3 x M 3 x [0, T] 

<C|x-y|(l + |ln|x-y||). (2.11) 

In particular the solutions of the ODEs (|2.7p and (|2.10p are uniquely defined 
as long as the distance between the plasma particles and the charges remains 
positive. 

Thanks to the hamiltonian structure of the system, the flow (X, V)(0, t) 
preserves Lebesgue's measure on M 3 xM 3 . As a result all the norms 
1 < p < +oo are preserved. Finally, it follows from (|2.5p . (|2.9p and from the 
fact that E is bounded that the density f(t) remains compactly supported 
for all t G [0,T]. 

We define the total energy associated to Problem (jl.ip by 

A' 



If 1 

H(f) = - \v\ 2 f{x, v) dxdv + - ^ \Va 

J a=l 



+ T,fr A ^dx+ 1 - ff P -^4dxdy+ 1 -y / -^—. 

\*~U *JJ \*~V\ 2£jlk-fcl 

One easily checks that for a solution to Problem (|1.1[) H(f(t)) is finite and 
constant on [0, T]. Due to the positivity of the interaction, the kinetic energy 
(|2.ip of the plasma part is bounded by H(f(t)) = H(f(0)). Therefore we 
may assume that the constant K appearing in (|2.4p does not depend on the 
time thanks to the conservation of the energy and of ||/||z,°°. 

Let us come back, for the moment, to the Vlasov-Poisson problem ig- 
noring the point charges. 
Denoting by 

P{t) = sup{M | (x,v) G aupp(/(t))} 
we conclude by (|2.4|) that 



P(t)<P(0)+K f P(s) 4 / 3 ds. (2.12) 
Jo 
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Indeed, the electric field computed on a characteristic X = X(t) is bounded 
by 



W x{t ^\ ± / \x(t)- x T dx = J \x{t)- x >v dx > 



hence ([ZH2D follows from fl27 

Obviously ()2.12p does not provide an a priori global bound for P(t). 
A refined estimate allowing to solve the 3-D problem has been obtained by 
considering the time averaging of the electric field along a trajectory, see [9j, 
[lUj . or [2J. The basic idea consists in partitioning the time interval [0, T] 
into small pieces (tj-i, U) for i = 1,2, . . . , of small length AT = |ij_x — 
For a given characteristic (X, V)(t), one writes thanks to Liouville's theorem 



1 1 dt\E(X(t),t) \<C I* dt f 



f(y,w, <i— i)t^77\ w7\i2 d 2/ ^ 



\X(t)-Y(t)\< 

(2.13) 

where (y,PF)(t) is a characteristic leaving at time U—x. 

Now, when 5^(t) is a trajectory of large velocity and large relative velocity 
at time (namely \w\ and \V(ti-\) — w\ are 0(P 4//3 )), we restrict our 
attention to the time integral 

dt 

(2.14) 



u _ x \X(t)-Y{t)r 

Here AT is chosen so small that the relative velocity remains large in that 
time interval {stability property). Then the time integral (|2.14p can be com- 
puted almost explicitely, using that X(t) — Y(t) essentially performs a free 
motion. As a consequence the contribution of (|2.14p to the time integral of 
the electric field is shown to be smaller than O(PAT) (see [9], [TO], [11] and 
[2] or Lemma [3] below). We call this contribution scattering plasma-plasma. 

The other contributions in (|2,13p can then be handled by means of static 
estimates relying essentially on Proposition [T] above. 

When a single point charge is present, the scenario changes dramatically, 
since the stability property for the trajectories of the plasma fails. Indeed 
the relative velocity of the plasma particles can change extremely fast if one 
of them collides with (or get very close to) the point charge. Nevertheless the 
interval of time in which the scattering charge-plasma takes place is so small 
that the contribution to the time integral (|2,14p can again be controlled (see 
Lemma H] below) . 

One also has to remark that, instead of using the maximal velocity Pit) 
as a control quantity, it is more convenient to use as in [T] the energy of 
a trajectory, the time derivative of which cancels the singular part of the 
electric field. 



6 



Once treated the single point charge-plasma in Section [31 we can turn 
to the iV-charges problem in Section 21 Indeed, the choice of AT and of all 
other parameters ensures that a plasma particle cannot get close to more 
than one point charge in every time interval As a consequence we 

can transfer the single charge analysis to the iV-charges problem with minor 
modifications. 



3 The plasma-charge model 

Let / = f(x,v,t) be the probability distribution of the plasma particles, 
and let £ and rj = £ denote the position and the velocity of a single point 
particle of unitary charge. Problem (|2.6p - (|2.1Up reads 



/ G L°° (L 1 n L°°) , p £ L°° (L 1 n L°°) , E = f^*p, (3.1) 



with 

f(X(x,v,0,t),V(x,v,0,t),t) = f(x,v,0), (t,i,»)eMxR 3 xI 3 , (3.2) 
where (X,V)(x,v,0,t) = (X,V)(t) satisfies 
'x(t) = V(t) 



\x(t)-m\ 

(X,V)(0) = (x,v), x^m 



V(t) = E (X(t),t) + ( 3 - 3 ) 



and 

U(t) = V (t) 

\v(t) = Em,t). 

The main result of this section is summarized in 



(3.4) 



Theorem 1. Let /o € L°° be a compactly supported probability distribution. 
Let (^Oj^o) £ R 3 x K 3 . Assume that there exists some 5q > such that 

mm{\x-£ \ | (x,v) € supp(/ )} > 5 . 

For all time T > there exists a unique solution (£, /) to Problem f)3 - 1 j) - 
(|3.4p on [0, T] with this initial datum. 



We shall prove that, assuming a solution to Problem (|3,ip - (|3.4p exists 
up to a fixed but arbitrary time T > 0, we have 

sup{\V(x,v,0,t)\ + — -3- —rr\ t€ [0, T] , (x, v) G supp(/ ) ) <C(T) 

[ \X{x,v,0,t) -£{t)\ J 

(3.5) 
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+ / T^Xdx+ 1 - ff P^EM dxd y. 



where C(T) is a constant depending only on /o and (Co,Vo)- As a conse- 
quence of (|3.5p the fact that a unique solution to Problem (|3.ip - (|3.4p does 
exist follows by rather standard arguments presented in Paragraph 13.31 at 
the end of the section. 

In what follows, Cj and K$,i = 0, 1, ... , will denote positive constants 
depending only on ||/o||l°° and H, where H denotes the global energy. 

In the case of one single point charge, the energy reduces to 
H =\J \v\ 2 f(x,v,t)dxdv + ^\ V (t)\ 2 

\x-£(t)\ °"~ ' 2JJ \x-y 
In particular, we have 

\ V (t)\ < V2H. (3.6) 
Following [1] we also introduce the pointwise energy of a plasma particle 

h(x,v,t) = b ~^ )|2 + 1 + K U (3.7) 

where K\ is a large constant. A possible choice is 

Kx > max(8#, 1). 

In particular, in view of (|3.6p this choice ensures that for all (x,v,t) 

|V(aj,u,0,t)| < 2Vh(X(x,v,0,t),V{x,v,0,t),t) . (3.8) 

As already mentionned, the energy turns out to be a relevant quantity to 
control the motion, since it controls both the velocity and the distance from 
£ of the characteristic under consideration. We remark that h is uniformly 
bounded on supp(/o) at time 0. 

In the following we shall use the short-hand notation (X(t),V(t)) = 
(X(x,v,0,t),V(x,v,0,t)) when the initial condition (x,v) is clear from the 
context. 

Differentiating along the characteristics of the plasma particles and using 
D-daaD we find 



h (X(t), V(t),t) = (V(t) - v (t)) ■ (E(X(t),t) - t)) 
from which 

j t Vh(X(t),V(t),t) | < \E(ffl,t)\ + \E(X(t),t)\. (3.9) 
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Note that the variation of h is controlled by the smooth part of the electric 
field and this is, of course, crucial. We introduce the quantity 

Q = Q T = sup {Vh(X(t),V(t) t t) | te [0,T],(x,v) Gsupp(/ Q )}. 

The remainder of this section is devoted to the proof of the following estimate 
from which (|3.5p follows immediately. 



Proposition 2. We have 

Qt< (Qo + CJexpiCitl + T)) VT>0. 

As explained earlier, the method relies on an suitable splitting of [0, T] 
into small intervals. More precisely, we set 

AT ' 



K 2 Qt 

where K 2 denotes a suitable large constant satisfying 

8K 1 

K 2 > 16 and — < -, 
K 2 a 

where K (depending only on ||/o||oo and H) is the constant appearing in 
Proposition [TJ 

Next, if AT < T we set 



n 



AT 



to = 0, t n = T, ti = iAT for i = 0, . . . , n — 1, 



so that 

n 

[0,T\ = \J[ti-i,ti] with |ij — ti-i\ < AT. 



i=l 

For i = 1 , . . . , n we define 



Qi =sup{v^(X(i),y(t),t) I t€ (ti_i,ti),(x,«) €supp(/(ti_i))} 

where here (X(t),V r (t)) = u, V(x, u, ii_i,t)) are the trajecto- 

ries at time t > U-i, leaving (x,v) at time £j_i. Finally, we set 

Qo = sup j\//i(x,7J,0) | (x, v) G supp(/ )| . 

In order to show Proposition [2] we will first establish the following basic 
inequality 

Proposition 3. Let T > such that AT < T. We have 
Qi < Qi-i + C 2 Q T AT, i = l,...,n. 



We claim that Proposition [2] follows from Proposition [3l Indeed, let us 
set To = 1/(AC 2 ). There are two cases. 

If ATq = 1/(K 2 Qt ) < Tq then Proposition [3] for To implies that for all 
i = 1, . . . , n 

Qi<Qo + C 2 iQ To AT <Q + 2C 2 T Q To , 

hence 

Qt < ° „ = 2Q . 

1 — 

Otherwise we have ATo > To, which means that 

Qt <^ = ±C 2 K 2 \ 

In both cases we obtain 

Qt <2Q + AC 2 K 2 ~ 1 , 

thus Proposition [5] holds up to time To. Let now T > Tq and k = [^]. 
Since Tq depends only on conserved quantities, we can iterate the previous 
arguments k + 1 times to get 

k 

Qt < Q( k +i)T < 2 h+1 Q + 4C 2 K^ 2? < 2 T ^ +1 (Q + 2C 2 K 2 l ) 

3=0 

and the conclusion follows. 

We now come to the main ingredients for proving Proposition We 
observe preliminary that, without loss of generality, we may assume 

Qi > K 3 > 1 (3.10) 

where K% is a constant depending only on K\ and K (therefore only on 
|| /|| 00 an d H) which will be specified in the course of the proof of Lemma 2] 
below. Indeed, otherwise we have 

Q l < K 3 < K 3 K 2 ^- < Qi_! + C 2 QAT 

provided that C 2 > K 3 K 2 , and Proposition [3] follows. 

Next, we notice that, by virtue of (|3.9p . (j3.8j) and (j2.4f) we have 

Vh (X(t), V(t),t) < Vh(x, v, ti-x) + 8KQ 4 / 3 AT (3.11) 

for t G and (x,v) 6 supp(/(tj_i)). Now, consider a trajectory 

(X(t),V(t)) satisfying 



Vh{X(t),V(t),t) =Q t 
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for some t 6 We then have by definition of AT 



Vh(x,v,ti~i) >Qi- -r^Q 1 / 3 > TT- 



(3.12) 



Therefore to control Qi it suffices to control the energy of those trajectories 
for which (l3T2|) holds. 

In order to bound the time integral on the right-hand side of (|3.9p we 
have to evaluate the integrals 



for high energy trajectories X(t), Y(t) which could possibly make very small 
the denominators in (|3.13p . There are various situations: 

(i) Both X and Y are far from £ (Lemma [3]), 

(ii) X is close to £ (Lemmas H] and [5]) , 

(iii) X is far from £ but Y is close to £, hence X and Y are far from each 
other (Lemma [7]). 

We shall handle each of these situations separately, achieving thereby the 
dynamical part of the proof. The remainder of the proof relies on rather 
straightforward estimates in phase-space. 

3.1 Preliminary estimates 

We start by establishing a lemma concerning the plasma-charge scattering. 
Lemma 1. For any (y,w) £ supp(/(ij_i)) we have, with (Y,W)(t) = 




and 




(3.13) 



(Y,W)(y,w,h-i,t) 




Proof. Setting t(t) 



\Y(t) 



£(i)|, we differentiate 




(3.14) 



then 




l(Y - Q • (w - n)} 2 

|Y-£| 3 



Proposition Q] and Cauchy-Schwarz inequality yield 



£>-^-8KQt 
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Therefore 



8K „ 1/3 



<\W- n\{ti) + \w- v\(ti-i) + -r^Q 

By definition of if 2 and since the first two terms in the right-hand side of 
the previous inequality are bounded by V2Qi we conclude the proof. □ 

We now introduce the quantities 

Ri = Q 3 / 4 and 5i = Q~ 7/8 . (3.15) 

Note that Ri corresponds to the maximal radius in the velocity space for 
which (|2.4p actually works, namely yields a linear estimate in Qi. 

As already mentionned, the choice of the parameters ensures stability 
for the quantity y/h. 

Lemma 2. Let (y, to) G supp(/(tj_i)) and (Y, W)(t) = (Y, W) (y, to, tj_i,t). 
If Vh(y,w,ti-i) > Ri then 

Vh (Y(t), W(t),t) > -Ri, Vt G [U-i,U]- ( 3 - 16 ) 

If Vh(y,w,ti-i) < Ri then 

Vh(Y(t),W(t),t)<^Ri, Vte^-i.ti]. (3.17) 

Proof. This is a straightforward consequence of (|3.9|) . observing that by 
choice of AT we have 8KQ* /3 AT < R4/2. □ 

The quantity Si is the radius of a protection sphere around the charge £, 
outside which the singular field created by ^ is relatively moderate, namely 

The next lemma deals with the plasma-plasma scattering when the in- 
fluence of the charge is not too large. The control of the corresponding time 
integral is the same as in [9], [TO], [11] and [2]. We shortly repeat it for 
completeness and also because, for the present problem, we need a different 
choice of the parameters. 

Lemma 3. Let li > 0. Assume that there exists a time interval 

J = (ti-i,ti) C (tj_i,tj) 

such that, for all t G J we have, with (X, V) (t) = (X,V) (x,v,ti—i,t) and 
(Y, W) (t) = (Y, W) (y, to, ij-i, i), where (x,v) and (y,w) G supp(/(fj_i)) 

mfmm{\X(t)-at)l\Y(t)-m\}>ti (3-18) 
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and 

\VQ i - 1 )-W(t i -i)\>Ri- (3-19) 

Then 

f U dt x{\X{t)-Y(t)\>k) 



Remark 1. We shall use this estimate with the choice Zj = Q, 



-2 



Proof. Let to £ \ii-i,U] be the minimizer of \X(t) — Y(t)\. Suppose for 
the moment that to £ (ti-iiU)- Setting D(t) = X(t) — Y(t), we have for 
t G (to,**) 

D(t) = D(t ) + D(t )(t - to) + p(t), (3.21) 

where 

p(t) = I ds(t-s)( 



- xjs) - m Y( S )-as) 

v,X( S )-e( S )| 3 \Y(s)-m\ 3 (3.22) 
+ E(X(s),s)-E(Y(s),s) 

By virtue of (|3.18p , Proposition [1] and the definition (|3.15p of R{ and 5, we 

have 

Ki)|<^-*o) (2Q] /4 + 8KQf 3 ) AT 

1 /. . v n 3/4 



hence 



< o(*-*o)<9- 



|/£>(t)|<^(t-to)fli. (3-23) 



On the other hand, by (IBTTHD and ([3~T9]) it holds 

\D(t )\ > \D(ti-i)\ - P \D(s)\ds 
>Ri- (2Q 7 / 4 + 8KQ 4 / 3 )AT 



hence 

\D(t )\ > ^. (3.24) 

We remark that the parameter 5i has precisely been chosen so as to 
ensure the above stability property, and that we have used the fact that Ki 
is sufficiently large with respect to K. 
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We have 

\D(t)\ > \D(t ) + D(t )(t-t )\-\p(t)\. 
Since -D(to) and D(to) are orthogonal, 

\D(to) + D(t )(t - t )\ 2 > \D(t )\ 2 + \D(t )\ 2 (t - t ) 2 > \D(t )\ 2 (t - t ) 2 . 

(3.25) 

Hence we have by ([3331) and (I3T241 

\D(to) + D(t )(t - to)\ > \D(to)\(t - t ) > -Riit - to) > 4|p(t)|, 
so that 

\D(t)\ 2 > ^\D(t )\ 2 (t - t ) 2 > ^R 2 (t - to) 2 . 



Finally, 



^ ^ X (\D(t)\>k) < r * x(\D(t)\ > k 



dt f 64 dt 



< 



3R 

a 



t0 <t<t +£t l 1 Jt>t +*t VR 2 (t - to) 2 



l%Ri 

For the integral on the time interval (tj_i,to) we use the time reversal. 
When to = M-i we use the same method, observing that D(to)-D(to) > 
so that (|3.25p still holds. Finally if to = U we use the time reversal. □ 

Next we need to control the charge-plasma scattering. Basically our 
aim is to show that the time spent by a trajectory in the protection sphere 
B(£(t), 5{) is very small. For proving this we apply the virial theorem argu- 
ment, introducing 

I(t) = \\Y(t)-m\ 2 - (3-26) 

Differentiating we get 

i = (Y-0-(W- V ) (3.27) 

and 

I = \ W — rj\ 2 + — 1- + (Y - ■ (E(Y) - E(0) ■ (3.28) 

Lemma 4. For (y,w) £ supp(/(tj_i)) ; suppose that 

Vh(y,w,ti-i) > B4. 
Consider (Y,W)(t) = (Y,W)(y,w,t i _ 1 ,t). Then the set 

J={te(U-i,u)\ \Y(t)-£{t)\<5i} 
is connected. Moreover, 

meas(J) < C 5 Q,~ 13/8 . (3.29) 
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Proof. Let to G J be a minimizer for I(t). By (|3.27p . 

\i(t)\<y/2I®\W(t)-7,{t)\, 

hence 

\Vi(t)\ < Q<. 

For t G [to,ti) we obtain 

y/w < vWoi + QjAt < ^Qr r/8 + 

therefore 

|/(t)| < 1. (3.30) 
Moreover, by Lemma [2] we have for t G [to,tj) 

v^(y(*),w(t),t)>^. 

Then by (^28^ and (EOl we have for t G [to, *i) 

/(t) > /i (Y(t), W(t),t) — K\ — \Y(t) - ttt)\\E(Y(t),t) - E(£(t),t)\ 
>±Rl- Kl -8KV2Qf 

= - A Qf-K 1 -8KV2Qf. 

Now, we observe that for K% > sufficiently large (depending only on K\ 
and K) we have by assumption (|3.10|) 

I{t) > l -P3 t . (3.31) 

Consider now (t_,t + ) C J a maximal connected component containing 
to- If to £ [t-,t+), /(to) > (if to = t+ we use the same argument via the 
time reversal). Then 

i(t)>i(t ) + -R?(t-t )>o, vtG [to,u). 

o 

Since I is increasing from to U P to ij, the trajectories cannot reenter in the 
protection sphere once escaped. Therefore J = (t_,t + ) is connected. 
Next, integrating twice (|3.3ip in time and using that /(to) > we get 

\% > I(t) > /(to) + ^Rf(t - to) 2 , Vt g J, 

so that 

(t - to) 2 < 8Q~ 13/4 , Vt G J, (3.32) 

and (|3.29p is proved. 

□ 
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We finally obtain the following variant of Lemma 0] 
Lemma 5. For (y,w) £ supp(/(ij_i)), suppose that 

Consider (Y,W)(t) = (Y,W)(y,w,t i _ 1 ,t). Then the set 

J = {te(U- 1 ,u)\\Y(t)-^t)\<2d i } 

is connected. Moreover, 

meas(J) < C 6 Q^ 15/8 . (3.33) 

Proof. It suffices to follow the proof of Lemma[Hstep by step, observing that 
estimate (|3.32p can even be improved if h(y,w,ti^\) > Qf/4. Indeed, we 
have in this case I(t) > Qf/8 for t E and everything goes exactly 

as before leading to (|3.33p . □ 

This concludes the dynamical preparation. We now come to the proof 
of Proposition [3l 

3.2 Proof of Proposition [3] 

In view of (|3.9p . in order to show Proposition [3] we need to control the 
integrals 

Ji= dt\E(£(t),t)\ and J 2 = dt\E (X(t),t) \. 
Lemma 6. We have 

Ji < C 7 QAT. 



Proof. Setting (Y,W)(t) = (Y,W)(y,w,t i ^ 1 ,t) we have 

Ji < dt dy dwf (y,w, U-i] 

Jti-i J (v,w)\\/h(y,ii),ti-i)<Ri 



1 



U-i J(y,w)\Vh(y,u>,ti-i)<Ri " " ^' ±7 — C(*)| 2 
+ / dydwf{y,w,ti-i) / dt — — — ^. 

Note that, by virtue of the stability property for yh (see Lemma [2]) and 
by ([3SD 

\W(t)\<3Ri, Vt€(ti_i,ti) 

when y/h(y,w,ti-i) < R{. Therefore, thanks to Liouville's theorem the first 
term in the right-hand side of (|3.34p can be controlled by 

[ dt I dydw f(y,w,t)— — \t-tk. 

Jti-i J\w\<3Ri \y -£(*)! 
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By virtue of Proposition [I] with R replaced by 3Ri this latter is bounded by 

C 8 Rt /3 AT = C s Q l AT. (3.35) 

We now turn to the second integral in (|3.34p . By virtue of Lemma [T] we 
can bound it by 

5Qi / dydw f(y,w,ti-i) 

J(y,w)\ h(y,w,ti-i)>Ri 

<5QiQJ 3/2 / dydwh(y,w,ti-i)f(y,w,ti-i) 



<c 9 q; 1/2 (h + i) 

< CxoQAT, 

where we have used that Qi > K3 in the last inequality. The conclusion 
follows. □ 

Lemma 7. Let (x,v) £ supp(/(ij_i)) and (X, V)(t) = (X,V)(x,v,ti-\,t) 
be a trajectory such that h(x,v,ti-i) > Then 

J 2 < CuQAT. 

Proof. Let (i_,i+) be the connected interval (if any) in which we have 
\X(t) — £(t)\ < 25{. By virtue of Proposition [T] and Lemma (estimate 

J dt \E (X(t),t) I < 4KQ 4 / 3 C 6 QT 15/8 < C 12 QAT. 



It remains to control the integrals 

r-t- fti 

'ti-i 



f dt\E(X(t),t)\ and f ' dt\E (X(t),t) |, 



which can be handled (using again the time reversal) in the same way. We 
write 



/ dt dydwf(y,w,ti-i) 
4 f 

V] / dtdydw f(y,w,ti 
~1 Js, 



X{t)-Y(t)\* 

1 



i=i " s > 



\X(t)-Y(t)\* 
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where 

51 = {{y,w,t)\ Vh(y,w,ti-i) < Ri}, 

5 2 = {(y,w,t)\ Vh(y,w,ti-i) > Ri and \X(t) - Y(t)\ < k}, 

5 3 = {(y,w,t)\ Vh(y,w,ti-i) > Ri, \X(t)-Y(t)\>k and \Y{t) - £(i)| < **} 
Si = {(y,w,t)\ Vh(y,w,t i _ 1 )> R u \X{t) -Y{t)\ > k and \Y{t) - £(i)| > 

and where 

k = Qf- 

Using stability for \fh (see Lemma [2]) as well as static estimates, the first 
integral J\ can be estimated as before in Lemma [6] (see (|3.35p ). so that 

Ji < C 13 QiAT. 

For the integral on S 2 , following [9], [10], [H] or [2] we get 

J 2 < t dt [ dy ._ P ^,. 2 < C u Q 3 ikAT < CuQiAT. 

Jti-i J\x(t)-y\<h \y-x{t)\ 



Next, for the integral ^3 we have 



J 3 < dydwf(y,w,ti-i) / dt 

Jh(v,w,U-i )>R? Jti-i 



X (\Y(t)-m\<Si) 



1 Hy ^U-,)>Rf a ~~ Jyai "4_r \Y(t)-X[t)\* 
Since \X(t) — £(i)| > 25i in (ij_i,t_) we have in S3 

\Y(t) - x(t)\ > \x(t) - m - \y(t) - m\ > ^. 

Hence, in view of the conservation of the energy and of Lemma H] (estimate 
<K29h ) we have 



1 



J3<j2 dydwf(y,w,U-i) dt X {\Y{t) - £(t)\ < 6t) 

°i Jh(y,«;,t i _ 1 )>R2 J ti _ 1 

< C 16 QAT. 

We finally estimate the last integral J^. By virtue of Lemma U for each 
(y, w) such that yfh(y,w } ti^\) > Ri we may split the set 

\Y{t)-m\>&i} 

into two at most intervals J l (y : w) and J 2 (y,w) for which 
inf \Y{t) -£(t)\>6i, k = 1,2. 

t£j k (y,w) 
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Note that at least one of the extremal points of J 1 or J 2 has to coincide 
with at least one of the £j_i or tj. Hence we have 



k=x Jh{y,w,U-i)>Ri 



dy dw f(y, w, t 



i-lj 



dt 



X (\X(t)-Y(t)\>h) 
\Y(t)-X(t)\* ■ 



It suffices to control the integral on J 1 (y, w) because the integral on J 2 (y, w) 
can be handled in the same way. 

We set J 1 (y, w) = (f_, f + ). Then we further split the integration domain 
as follows 



{(y,w) | >/%,«;, ti_i) > Ri} = U Sf\ 



where 



SF = {(y,w)\ Vh(y,w,ti-i) > Ri and \W(t-) -V(t-)\ < Ri}, 
{(y,w)\ Vh(y,w,t i - 1 )>R i and \W(t-) - V(t„)\ > R4}. 



:(2) 



We recall that W{t) denotes the velocity of the plasma particle leaving (y, w) 
at time 

First, for (y,w) G S4 we have 



\W(t)-V(t)\<-R h Vt€Ji(y,w). 

Indeed, both X(t) and Y(t) remain at least at distance Si of the charge in 
Ji(y,w), therefore 



<20 7/4 +8A'«* /3 , 



so our choice of AT ensures stability for the velocities (see the proof of 
Lemma E]). Hence, applying Liouville's theorem we obtain 



(i) 



dydwf(y,w,t 



i-l 



dt 



X (\X(t)-Y(t)\>h) 



M \Y(t)-X(t)\ 
< I dydwf(y,w,ti_i) 

dt I 

\w-V(t)\<3Ri/2 



U dt xi\V(t)-W(t)\<3R t /2) 



< 



dt 

i-l 

< C 17 QiAT. 



U—i 

dydwf(y,w, t) 



\Y(t)-X(t)\* 
1 



\y-x(t)\* 



where we have used Proposition [TJ in the last inequality. 
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Finally, let (y,w) G sf \ Since \X(t) - f(t)| > 2<5j > Si on J l (y, w), we 
may now apply Lemma to obtain 

r X {\x{t)-Y{t)\>k) ^ Ci 8 

Since = Q~ 2 this yields 

This concludes the proof of the Lemma. □ 

Combining Lemmas [6] and [7] we may finally turn to the 
Proof of Proposition [3] completed. 

Let (X,V)(t) = (X, V) (x,v,ti-\,t) be a trajectory such that 
Vh(X(t),V(t),t) = Qi for some te[ti-i,ti]. 
By virtue of ()3.12p we have 

y/h(x,v,ti-i) > 

On the other hand, integrating first (|3.9p on and applying then 

Lemmas and to the high-energy trajectory (X, V) we obtain 

Vh(X(t),V(t),t) < Vh(x,v,t i - 1 )+ t dt (\E(X{t),t)\ + \E(£(t),t)\) 

Jti-l 

< Qi-i + C 20 QAT, 

whence 

Q l <Q l -i + C 20 QAT 

and the conclusion follows. 

3.3 Proof of Theorem Q] completed 

We finally complete the proof of Theorem Q] in the following way. We first 
regularize the Coulomb kernel l/\x\ at a small distance from the origin, say 
e, obtaining a solution (£ e (t), rj £ (t), f £ (t), X £ (t), V £ (t)) of the corresponding 
regularized version of Problem (|3.ip - (|3.4|) . The fact that the corresponding 
global energy H £ is uniformly bounded by H provides uniform bounds for 
the kinetic energy (|2.ip of f £ and for \r/ £ \ (see (|3.6p ). 
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We choose e smaller than 1/C(T), where C(T) is the a priori bound 
appearing in (|3.5p . Next, we take T £ maximal such that 

min{|X £ (x,?;,0,t) -rf(t)\ | G supp(/o)} > £, i€[0,T e ). 

On [0,T e ), the extra field created by the charge coincides with the one 
of Problem (jl.ip for all trajectory starting from supp(/o). Therefore the 
previous analysis applies, yielding (|3.5|) , We conclude that T £ = T, so that 
(|3.5p provides uniform L°° bounds for p e and E £ on [0, T\. It follows that 
(^ e , n £ ) is uniformly bounded and equicontinuous on [0, T] and that {X s , V s ) 
is uniformly bounded and equicontinuous on supp(/o) x [0,T]. Hence one 
easily passes to the limit e — > to get existence of a global solution (£, /) 
satisfying the desired conditions. 

Finally, uniqueness is achieved by means of a Gronwall inequality, using 
almost-Lipschitz regularity for E (see (|2.11|) ) and the lower bound (|2.9p 
together with the assumption on the support of /q. We omit the details. 

4 The case of many point charges 

The purpose of the present section is to extend the existence and uniqueness 
result of the previous section to the case of many point charges. We will 
establish the following 

Theorem 2. Let N > 1. Let fo £ L°° be a compactly supported distribution 
and let {£ q0 , ??ao}a=i such that £ q0 ¥= £/30 fora^/3. 
Assume that there exists some 5o > such that 

min{|x - £qo| | (x,v) £ supp(/ )} > S . 

Then for allT > 0, there exists a unique solution to Problem (|2.6|) - (|2.10D 
on [0, T] with these initial data. 

We proceed as in the proof of Theorem [TJ considering a solution on [0, T] 
obtained by regularizing the singular interaction kernel in Problem (jl.ip and 
establishing estimates depending only on H/olloo and H. Here again, C, Cj 
or Ki, i = 1, . . . , denote positive constants depending only these quantities, 
C changing possibly from a line to another. 

First, we infer from the conservation of the total energy H(f) that 

!&(*)-&(*)!> A, Vie[0,T], (4.1) 
where A = l/(2H). Also, 

\rj a (t)\<V2H, Vte[0,T]. (4.2) 
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For a = 1, . . . , N we introduce the energy related to the a-th charge 

h a {x,v,t) = \\v - ri a {t)\ 2 + —— + Ki, (4.3) 

2 \x-£ a {t)\ 

where K\ has already been defined in the previous section. 
We set 

Q= max sup-j \/h^(X{x, v,0, t),V(x,v, 0,t),t) \ t G [0, T], (x, v) G supp(/o) > ■ 

a=l,...,N I ' J 

As before, we split the interval [0, T] into small intervals of length 

smaller than AT, where 

AT = — . 

Here K 2 , depending only on H and ||/o||ooj is chosen in a similar way as in 
the previous section, and we assume moreover that 

K 2 > ^. 

A 

Therefore, if (X, V)(t) = (X, V)(x,v,ti_i,t) denotes a plasma particle leav- 
ing (x,v) at time we have thanks to (|3.14p 



j t \x(t)-ut)\ 



< V2Q < (4.4) 
- ^ - 48 AT v ; 



so that for all a = 1, . . . , N 



x € S ^) 41(^5^),^ Vt€[t<_i,ti] (4.5) 

and 

xeB^),^ ^I(i)GB^(t),^j ViG [ti-i,**]. (4.6) 

In view of (|4,ip , equation (|4.5p means that a plasma particle starting close to 
the a-th charge at time cannot approach the other charges on 
This property will enable us to isolate each charge and to apply the previous 
analysis to the present case. 

For i = 1, . . . , n we define 
Qi= max sup\ y^(X(t),V(t),t) \ t e (ti-i,ti),(x,v) e supp(f(ti^i))\ 

a=l,...,N K J 

where (X,V)(t) = (X, V) (x,v,ti-i,t). Finally we set 

Qo = max sup \ vh(x, v, 0) | (x, v) G supp(/o) \ 

a=l,...,N I J 



22 



so that 

Q = max Qi. 

i=0,...,n 

As explained in the previous section, Theorem [2] is a consequence of the 
following variant of Proposition [3] for the present case 

Proposition 4. Let T > such that AT < T. For all i = 1, . . . , n we have 

Qi <Qi-i + CiQAT. 

We will only present the proof of Proposition [H Again we may assume 
that 

Q l > K 3 > 1. 
We introduce the security balls 

B^bU^x),^), /3 = 1,...,N (4.7) 

and the complementary set 

N 

n = m 3 \ |j Bp. 

Next, we set for t G 

fp(x,v,t) = XB f(x,v,t) and f{x,v,t)=xnf{x,v,t), (4.8) 

so that 

N 
/3=1 

and we denote by Ep and E the corresponding electric fields, so that 

N 

E = J2 E P + E- 

(5=1 

Lemma 8. Let a € {1, ... , iV}. Xei (x, t>) G supp(/(ij_i)) such that 
x e B a and h a (x,v,ti-x) > -rQi- 

Then 

y/h^(X(t),V(t),t) < s/h^(x,v,t i - 1 ) + C 2 QAT, te [U-uh). 
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Proof. Introducing the field 



we compute 
d 



(// v h a (x(t),V(t),t) 



<\E(X(t),t)\ + \E(Ut),t)\ 

(4.9) 

+ \F a (X{t),t)\ + \F a (Ut),t)\- 



Thanks to (|4.5p and (|4.ip . we observe that X(t) remains far from the 
charges £p(t) for f3 ^ a. Hence the fields F a (^ a ) and F a (X) are bounded on 
[ii-l,tj] and we have 

1 ' dt (\F a (&(t),t) | + \F a (X(t),t) I) < CAT < CQAT. (4.10) 

Jti-i 

Moreover, since E + Y^p^a^fi + F a is bounded by 0(Q^ 3 ) away from 
the o-th charge, we may follow step by step the dynamical preparation 
performed in the previous section, simply replacing the quantities £(t), Vh(t) 
and J(t) by £ a (t) = \X(t) - £ a (t)\, VK(t) and I a (t) = \X{t) - &(t)| 2 /2. In 
particular, adapting the proofs of Lemmas [6] and [7] we find 

t dt(\E a (Ut),t)\ + \E a (X(t),t)\)<CQAT. (4.11) 

J ti-l 

It remains to estimate the contributions of the other fields. Let y G Bp. 
Thanks to (jiTT]) and (|i"3j) we have 

1 7A 

\X(t) - Y(t)\ > \Ut) ~ m\ - \X(t) - Ut)\ - \Y(t) - > — , 

hence Ep(X) is bounded on By (|4.ip and (|4.6p Ep{^ a ) is also 

bounded. Therefore 

^ tft (|^ (Ut),t) | + |£g (X(t),i) |) < CQAT. (4.12) 

ti " 1 

We finally estimate the contribution of E. By (|4.ip and (|4.6p E{^p) is 
bounded, thus 



/"' (E 1 (£ a (t),t) | < CAT < CQAT. (4.13) 



In order to estimate E(X(t),t) we distinguish between two cases. 
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We assume first that x G B(£ a (t{-i), A/16). For t/ 6 we then have 
\X(t) - Y(t)\ > A/48 by flUD and ([4~4]h Hence £(X) is bounded and we 
obtain 

£ dt \E (X(t),t) | < CQAT if x € B ^) • (4-14) 

Otherwise, we have x G i? a \ i?(£ a (tj_i), A/16). Let y £ £1. In view of 
gZEJ the particles X(t) and F(t) both remain in R 3 \ UpB^t), A/24) on 
[ti-i,ti]. One may therefore neglect the plasma-charge interaction and only 
take into account the plasma-plasma interaction. Following the arguments 
in the case without charges (see [H]), or adapting Lemmas [3] and [7] we obtain 

jf' dt\E(X(t),t)\<CQAT if x eB a \B(^ a (ti-i),^y (4.15) 

Gathering estimates (|4.10p to (|4.15f) and using (|4.9p we are led to the 
conclusion of Proposition □ 

Our next result concerns the variation of the energies h a away from the 
charges. 

Lemma 9. Let (x,v) G supp(/(tj_i)) such that x G £1. 
For a// a G {1, . . . , N} we have 



h a {X(t),V(t),t) < ^h a {x,v,t^x) + C 3 QAT, t G (ti_i,ii). 

Proof. We have Z(t) G K 3 \ U a 5(^ a (t), 5A/48). Imitating the proof of 
Lemma El it only remains to estimate 

U 

dt\E a (X(t),t) | 

ti-l 

for all a. We proceed as before, writing 

B a =B Ua{ti-l), ^\ U ^B a \ B Ua{ti-l), ^ 

noticing that \X(t) — Y(t)\ > A/48 when y belongs to the first set and 
ignoring the effect of the charges in the second set. The conclusion follows. 

□ 

Finally, we can compare the energies as follows 

Lemma 10. There exists a constant C4 such that for t G [0, T] we have for 
all a, j3 



hf,(X,V,t)<y/h^(X,V,t)+Ct, VX G B U a (t), VV G 
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Proof. Our aim is to show that 
2 , ,o .. 2 



+ |^| 2 - 2V • ^ < — - + |t/ q | 2 - 2F • r? Q + 2C| + V2/^C 4 , 



which holds if 



\v\< ^- 1 - i (\v a \ 2 -\vM 2 + 



2\ria-Vp\ V P l^-Cal l^-Cjsl 



+ 2C| + 72^C 4 ^ . 



We have |X - ^| > \X - £ a \ by virtue of flU}. Using ([32]) and we 
find that the right-hand side in the above inequality is larger than 

' ' -2H + 2Cj + 2\V\C 4 ) , 



4V2i7 

which is larger than \V\ provided that C 4 > 2\f2M. □ 

We may finally turn to the 
Proof of Proposition [4] 

Let t G [ti-i, ij] and /3 G {1, . . . , N} such that 

Qi = hp {X(t),V(t),t) ■ 
There are three possibilities. 

If x G Bp then X(t) G B(£p(t),7\/48) on [ii-;L,ii]. In view of gJl, sta- 
bility for holds on [£j_i, U]. In particular, we may choose K2 sufficiently 
large so that y/hp{x, v, ij-i) > Qj/2. Hence Lemma [8] yields 

Qi<Qi-i + C 2 QAT. (4.16) 

If x G we use Lemma [9] for vfya to obtain 

< Qi-i + C 3 QAT. (4.17) 

Finally, if x G -B Q for some a ^ (3 then stability for v^a holds on [ti_i, t»] 
by ()4.9p . Since by Lemma [10] at time t we have 

7^(X(t), V(t),t) > y^(X(t), V(t),t) ~C 4 >^ 

provided that Qi > K3 > AC 4, we may choose -K2 such that at time tj_i 
it holds y/h^(x, v, U-i) > Qi/2. Hence we may apply Lemma [8] to \/7ta- 
Relying then once more on Lemma [TU] at time t we get 

Qi < Qi-i + C 4 + C 2 QAT = Qi_! + (C 4 ^2 + C 2 )QAT. (4.18) 

Setting Ci = max(C2, C3, C4K2 + C2) the conclusion follows from (|4.16p . 
dHHH and (|Qg]> . 
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5 Final remarks and comments 



The first comment concerns the regularity of the solution we have con- 
structed in Theorems [1] and [2j The fact that t i— > f(t) propagates the C k 
regularity of the initial condition is standard; this is a consequence of the 
almost-Lipschitz regularity of the electric field E (see (|2.1ip ). In particular 
if fo £ C 1 the solution we have constructed is a classical solution to Problem 

(HP- 

The second comment regards our hypotheses. We assumed the charges 
with an island of vacuum and /o compactly supported in v . Of course the 
total energy can be finite also in absence of these hypotheses. However 
when the trajectories can be arbitrarily close to the charges or when they 
have arbitrarily large velocity, there is not any uniform upper bound for the 
initial energy of all the trajectories Qq. Therefore our method requires new 
technical ideas to take into account these large tails. 

Finally we mention that the case of an interaction plasma-charge of dif- 
ferent signs eludes our techniques: new ideas and estimates are needed. 
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